Pulsed-gate measurements of the singlet-triplet relaxation time in a two-electron 

double quantum dot 
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A pulsed-gate technique with charge sensing is used to measure the singlet-triplet relaxation 
time for nearly-degenerate spin states in a two-electron double quantum dot. Transitions from the 
(1,1) charge occupancy state to the (0,2) state, measured as a function of pulse cycle duration and 
magnetic field, allow the (1,1) singlet-triplet relaxation time (>70 /xs) and the (0,2) singlet-triplet 
splitting to be measured. The use of charge sensing rather than current measurement allows long 
relaxation times to be readily probed. 

PACS numbers: 73.21. La, 73.23.Hk, 85.35.Gv 



Semiconductor quantum dots are promising systems 
for the manipulation of electron spin because of the rela- 
tive ease of confining and measuring single electrons [J . 
In order to make use of the spin degree of freedom as a 
holder of either classical or quantum information, it is 
first necessary to understand and characterize the mech- 
anisms that lead to spin relaxation and decoherence. 

Previous studies of spin relaxation in quantum dots 
have focused on systems with large energy splittings of 
the relevant spin states, either singlet and triplet two- 
electron states, split by ~ 600 /icV 2], or Zeeman states, 
split by >200 /xeV Ulli]. In these cases, spin-orbit 
mediated processes are expected to dominate spin re- 
laxation '0, '?] , while processes involving hyperfine cou- 
pling of electron spin to nuclei are relatively ineffective for 
splittings greater than the hyperfine coupling scale ~0.1 
meV /^/ITe^I MgV, where Ne is the effective number of 
nuclei interacting with the electron spin 9, 10] . Spin 
relaxation between nearly degenerate electronic states is 
particularly relevant to the problem of controlled entan- 
glement, as singlet-triplet splitting goes to zero as the 
entangled spins become spatially separated. 

In this Letter, we present pulsed-gate measurements 
of the singlet-triplet relaxation times and splittings in 
a two-electron double quantum dot. Spin states of the 
weakly-coupled (1,1) charge configuration (integer pairs 
specify the equilibrium charge occupancy on the left and 
right dot) are nearly degenerate. A spin selection rule im- 
pedes the transition from (1,1) spin triplets to the (0,2) 
singlet. To probe spin relaxation, the double dot is ini- 
tialized, then loaded into the (1,1) state, then placed in 
a configuration where (0,2) is the preferred charge dis- 
tribution in equilibrium. This transition probability is 
measured from a calibrated charge sensing signal as a 
function of the period of this three-step cycle. The result- 
ing characteristic time of ^70 /iS sets the time scale for of 
the spin-blockaded relaxation channels from (1,1) to the 
singlet state of (0,2). Measurements of the (1,1) to (0,2) 
transition probability as a function of detuning, e, show a 
strong dependence on perpendicular magnetic field, B±, 



reflecting the singlet-triplet splitting of the (0,2) charge 
state, J as discussed below. 

Samples are fabricated from a GaAs/Alg.sGao.rAs het- 
erostructure grown by molecular beam epitaxy. Electron 
beam lithography and liftoff are used to create Ti/Au 
gates that deplete a 100 nm deep two-dimensional elec- 
tron gas with electron density 2xl0^^cm^^ and mobility 
2x lO^cm^/V-s. Gates 2-6 and 12 form the double quan- 
tum dot [see Fig. 1(a)]. Gates 2 and 6 are connected via 
bias tees to dc voltage sources and to pulse generators 
through coax cables with ~20 dB of inline attenuation 
[Til ]. A QPC charge detector is created by depleting gate 
1. Gates 7-11 are unused. The QPC conductance, Gs, 
is measured using standard lock-in amplifier techniques 



a) 9^ 1011 12 



b) -5 dGs/dV, (a.u.) 





^50 V2(mV) -400 




-415 V2(mV) -400 



-415 Vj (mV) -400 



FIG. 1: (a) SEM image of a device identical in design to the 
one used in this experiment. Gates 2-6 and 12 define the 
double dot. • denotes an ohmic contact, (b) Large-scale plot 
of dGs/dVe as a function of V2 and VJj. Charge states are 
labelled {M,N), where M{N) is the time averaged number 
of electrons on the left (right) dot. (c) Zoom in near the 
(1,1) to (0,2) charge transition, (d) 25 mV pulses with a 50% 
duty cycle and 10 pis period are applied to the coax connected 
to gate 2. This results in two copies of the charge stability 
diagram shifted relative to one another hy AV . 
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with a 1 nA current bias at 93 Hz. The electron tem- 
perature, Te^l35 niK, was determined from Coulomb 
blockade peak widths. We present results from a single 
sample; a second sample with slightly larger lithographic 
dimensions gave qualitatively similar results 12J. 

QPC charge sensing is used to determine the abso- 
lute number of electrons in the double dot. Figure 1(b) 
shows a large-scale charge stability diagram for the dou- 
ble dot. As electrons enter or leave the double dot, or 
transfer from one dot to the other, Gs changes, resulting 
in s harp features in dGs/dVe (numerically differentiated) 
[T^ll3| . In the lower left corner of Fig. 1(b), the double 
dot is completely empty. As the gate voltages are made 
more positive, electrons are added to the double dot. We 
will focus on the two-electron regime near the (1,1) to 
(0,2) charge transition [Fig. 1(c)]. In this location, charge 
transport in a double dot shows a striking asymmetry in 
bias voltage due to spin selection rules (Pauli blocking) 
[13. HH- At forward bias, transitions from the (0,2)s to 
the (1,1) singlet state (l,l)s are allowed. However, for 
reverse bias, the (1,1) to (0,2) transition can be blocked 
if the (1,1) state forms a (1,1)t because the (0,2)t state 
resides outside the transport window due to the large 
singlet-triplet splitting in (0,2). This asymmetry results 
in current rectification, which is used in present pulsed 
gate measurements. 

The double dot is electrostatically driven by apply- 
ing pulses to gates 2 and 6. Figure 1(d) shows a charge 
stability diagram acquired with square pulses applied to 
gate 2 (F/=25 mV, 50% duty cycle, period r=10 fis). 
This results in two copies of the charge stability diagram, 
the right-most (left-most) charge stability diagram re- 
flects the ground state charge configuration during the 
low (high) stage of the pulse sequence. The gate-voltage 
offset between the charge stability diagrams, AV, is used 
to calibrate pulse amplitudes. Additional calibrations are 
performed for gate 6, which primarily shifts the honey- 
comb in the vertical direction (not shown). A linear com- 
bination of pulses on gates 2 and 6 can be used to shift 
the stability diagram in any direction in gate space. 

As a control, we compare a forward pulse sequence, 
where spin selection rules are expected to be important, 
with a reversed pulse sequence that does not involve spin 
selective transitions. The forward pulse sequence (Fig. 
2(a)) begins with the gates at point E for 10% of the pe- 
riod, emptying the second electron from the double dot, 
leaving the (0,1) charge state. The gates then shift to 
the reset point R for the next 10% of the period, which 
initializes the system into the (1,1) configuration. The 
interdot tunnel coupling, t, is tuned with t<kBT so that 
the (1,1) singlet-triplet splitting j^At2/U<^kBT , where 
U is the single dot charging energy. Due to this degener- 
acy, and at low fields such that \g\^BB<kBT (|,g|^0.44 
for GaAs), we expect to load into (l,l)s or any of the 
three (1,1)t states with equal probability. For the fi- 
nal 80% of the period, the gates are at the measurement 
point M where (0,2)s is the ground state. Figure 2(b) 
illustrates the possible (1,1) to (0,2) transitions. If the R 
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FIG. 2: (a) Sensor conductance, Gs, as a function of V2 and 
Ve while applying the forward pulse sequence (see text) with 
T=10 /is and i3x = 100 mT. Spin-blocked transitions result in 
some (1,1) charge signal in the (0,2) pulse triangle (bounded 
by the red lines). Outside of the pulse triangle it is possible 
to access other charge states. This relaxes the spin blockade 
(see color-coded level diagrams), (b) Level diagram illustrat- 
ing the possible transitions from (1,1) to (0,2). Black arrows 
indicate fast transitions, grey arrows indicate spin blocked 
transitions, (c) Gs as a function of V2 and Ve while applying 
the reverse "control" pulse sequence with r=10 /iS and B± 
= 100 mT. The (0,2)s to (l,l)s transition is not spin blocked. 
As a result, there is no detectable pulse signal in the pulse 
triangle (bounded by the red lines), (d) Level diagram il- 
lustrating the (0,2) to (1,1) transition. A best-fit plane has 
been subtracted from the data in (a),(c) to remove signal from 
direct gate to QPC coupling. 



step loads the (l,l)s state, tunneling to (0,2)s occurs on a 
timescale given by the interdot tunneling rate, r(e) (We 
estimate the slowest r(e)(l /Lts)~^ from finite bias data 
01). If the ms=0 (1,1) triplet state (1,1)to is loaded, 
it dephases into (l,l),s on a timescale of T2 (expected 
to be <100 ns [II 111 113) followed by a direct transi- 
tion to (0,2)3. About half the time the R step will load 
the ms=l (1,1) triplet state (1,1)t+ or the ms=-l (1,1) 
triplet state (1,1)t-- At low S^, (0,2)t is inaccessible, 
and a transition from (1,1)t+ or (1,1)t- to (0,2) requires 
a spin flip and will be blocked for times shorter than the 
singlet-triplet relaxation time tst- 

In Fig. 2(a) the average charge sensor signal, Gs, is 
measured as a function of the dc gate voltages V2 and 
Vg, while a pulse sequence is repeated. This has the ef- 
fect of translating the points E, R, and M throughout the 
charge stability diagram, keeping their relative positions 
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constant. Because most of the time is spent at point M, 
the grayscale data primarily map out the ground state 
population for this point, with plateaus at Gs^O.O, 6.0, 
16, and 23 xlO~^e^/h indicating full population of the 
(1,2), (0,2), (1,1), and (0,1) charge states respectively. 
The pulse data differs from ground state data only when 
point M resides in the triangle defined by the (1,1) to 
(0,2) ground state transition and the extensions of the 
(1,1) to (0,1) and (1,1) to (1,2) ground state transitions 
(bounded by the red marks in Fig. 2(a)). Within this 
"pulse triangle" transitions from (1,1) to (0,2) may be 
blocked as described above, and the charge sensor reg- 
isters a conductance intermediate between the (1,1) and 
(0,2) plateaus. If M moves above the pulse triangle (red 
dot in Fig. 2(a)), the (1,1) to (0,2) transition can occur 
sequentially via (1,2) with no interdot tunneling: a new 
electron enters the right dot, then the electron in the left 
dot leaves. Likewise, if M moves below the pulse trian- 
gle (orange dot in Fig. 2(a)) the transition can occur via 
(0,1): the left-dot electron leaves, then a new electron 
enters the right dot. By similar logic, point R must be 
to the left of the (0,1) to (0,2) transition extension (dot- 
ted line in Fig. 2(a)) to avoid resetting through (0,2) and 
preferentially loading (l,l)s- Figure 2(a) shows a signal 
of llx lO^'^e^/h in the pulse triangle for r=10 /iS, which 
indicates that approximately 50% of the time the dots re- 
main in (1,1) even though (0,2) is the ground state. This 
is direct evidence of spin-blocked (1,1) to (0,2) transi- 
tions. 

If the the pulse sequence is reversed, so that the reset 
position R occurs in (0,2), where only the singlet state is 
accessible, and M occurs in (1,1), tunneling from R to M 
should always proceed on a time scale set by the interdot 
tunnel coupling, since the (0,2)s to (l,l)s transition is 
not spin blocked. As anticipated, no signal is seen in the 
pulse triangle for this reversed "control" sequence (Fig. 
2(c)). 

Spin selectivity of the forward pulse sequence in Fig. 3 
can be used to measure J as a function of B± 18] . Figure 
3(a) shows Gs as a function of V2 and Vq while applying 
the forward pulse sequence with B±=1.2 T and t=10 /iS. 
For these data, (0,2)t resides outside of the pulse triangle 
{JEm, the mutual charging energy) and the (1,1)t to 
(0,2) transitions are spin blocked. For B±=1A T [Fig. 
3(b)] the (0,2)t state is low enough in energy that the 
(1,1)t states can directly tunnel to the (0,2)t manifold at 
high detunings. Now (1,1) to (0,2) tunneling can proceed, 
and there is no longer a (1,1) charge signal in the (0,2) 
region of the pulse triangle at high detuning. This cuts 
off the tip of the pulse triangle. The spin-blocked region 
continues to shrink as B± is increased. From these data, 
we find J-340, 280, and 180 /zeV for B±=1A, 1.6, and 
1.8 T, respectively ^19]. 

The time dependence of the charge sensing signal can 
be investigated by varying r, the overall period of the 
cycle. Figure 4(a) shows Gs as a function of V2 and 
Ve acquired using the forward pulse sequence with r=8 
/iS at i?_L=100 mT. A clear pulse signal is observed in 
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FIG. 3: Forward pulse sequence results with r=10 fis. (a-d) 
Gs as a function of V2 and Ve for increasing B±. B± reduces 
J. Eventually (0,2)t is lowered into the "M" pulse window 
[inset of (b)]. At this point, (1,1)t to (0,2)t transitions are 
energetically possible and the transition from the (1,1) to (0,2) 
charge state is no longer "spin blocked" . This cuts off the tip 
of the pulse triangle in (0,2), see (b). A best-fit plane has 
been subtracted from the data in (a)-(d). 



the pulse triangle. As r is increased, the pulse signal 
decreases as shown in (b-c). Gs is measured inside 
the pulse triangle {V2,Via held fixed at -403,-523.8 mV, 
respectively) and is plotted as a function of r in Fig. 
4(d). In (1,1), Gs--20xl0"V/h, whereas outside the 
pulse triangle in (0,2), Gs'^lOx 10~^e^/h. For small r, 
G5'--^15x 10~^e^/h in the pulse triangle. At long r, Gs 
approaches lOxlO^'^e^/h in the pulse triangle, which in- 
dicates complete transfer from the (1,1) to (0,2) charge 
state. 

These data are consistent with spin-blocked transitions 
preventing the (1,1) to (0,2) charge transition. Approxi- 
mately 50% of the time, the (1,1) R pulse loads into ei- 
ther (1,1)t-i- or (1,1)t-. These states may relax into the 
(l,l)s state and then tunnel to (0,2)g on a timescale set 
by TsT- For T<^TST, we expect that 50% of the (1,1) to 
(0,2) transitions will be spin-blocked, resulting in ~50% 
(1,1) pulse signal in the (0,2) pulse triangle. For ttst, 
the (1,1)t-i- and (1,1)t- states have ample time to relax 
to (l,l)s, after which a (l,l)s to (0,2)s transition can 
take place. Thus, for r long compared to tst, the pulse 
signal to approaches the (0,2) level, indicating full trans- 
fer from (1,1) to (0,2). In the intermediate t regime, the 
sensor signal due to spin-blocked transitions decays as a 
function of time on a timescale that is characteristic of 

T ST- 

The experimental data in Fig. 4(d) are fit assuming 
exponential singlet-triplet relaxation. Due to the slow 
measurement rate of the charge sensor (~100 ms), Gs 
is proportional to the time-averaged occupation of the 
left dot. Modeling an exponential decay of the sensing 
signal weighted over the 80% of the cycle corresponding 
to the (0,2) measurement gives Gs{t)=A+B{tst/t){^ — 
g-o.8r/rsT^^ whcrc A is the conductance asymptote at 
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FIG. 4: Forward pulse sequence results with _Bx = 100 mT. 
{a)-(c) Gs as a function of V2 and Va for increasing pulse pe- 
riods, T. For longer periods, singlet-triplet relaxation occurs 
and the (1,1) to (0,2) transition proceeds, reducing the signal 
in the (0,2) pulse triangle, (d) Gs with l/2=-403.0 mV, Ve=- 
523.8 mV (gate voltage position indicated by the black trian- 
gle in (c)) measured as a function of r. Gs^lOx 10~^e'^/h in 
the (0,2) charge state, while G3=20x 10"^e^/h in the (1,1) 
charge state. For short r, the (1,1) to (0,2) transition is 
blocked approximately half of the time, resulting in a pulse 
signal of 15xl0"^e^/h in the (0,2) pulse triangle. The (1,1) to 
(0,2) transition probability increases (sensor signal decreases) 
with T due to spin-relaxation with a characteristic time scale 
of 70±10 /IS. A best-fit plane has been subtracted from the 
data in (a)-(c). 



(1,1) to (0,2) transition, tst decreases, to 20±5 /iS at 
V2=-403.8 mV and y6=-523.0 mV. Closer to the tip of 
the pulse triangle, tst decreases due to thermally acti- 
vated exchange with the leads (See Fig. 2, red and orange 
diagrams), thus the 70 /iS relaxation time represents a 
lower bound on the spin relaxation time within the (1,1) 
manifold 

In GaAs quantum dots, spin-orbit and hyperfine in- 
teractions are expected to be the dominant mechanisms 
for spin relaxation. The spin-orbit contribution to Ti 
depends strongly on B and theory predicts Ti=10 ms 
for B=l T 7]. The hyperfine spin relaxation rate is 

THF^rph(^^)2j^f^ where Tph-^xlO^-^ reflects 

the electron-phonon coupling, i?„=0.135 meV is the 
nuclear energy for a fully polarized system, Gc= 1.25 
is derived from the correlation of the nuclear spins, 
7n=0.045 depends on wavefunction overlap, the nuclear 
spin 1=3/2, and iVe^lO^ 9]. Our sample is tuned with 
j^4t2/C/<Cfcsr. In this case, the Zeeman splitting be- 
tween the TOs=0 state and TOs=±1 states sets the rele- 
vant energy scale. Using j^2.5fieV gives 1/Thf~^5 MS, 
roughly consistent with the experimental results. The 
relaxation times found here can be compared with those 
obtained by Fujisawa et al. who measured rsT=200 /iS 
in a vertical dot 0. In those measurements, the sizable 
singlet-triplet splitting implies that spin-orbit processes 
rather than nuclear processes dominant, in contrast to 
the present experiments. 



long times (full occupation of the (0,2) state) and B is 
the additional conductance in a short pulse due to the 
blocked states (approximately 50% of the (0,2) to (1,1) 
step height). The best fit to the data in Fig. 4(d) gives 
A=0.009 e^/h and B=0.007 e^/h, consistent with these 
expectations. In the center of the pulse triangle the best- 
fit TST reaches a maximum of tst=70±10 fis. Near the 
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